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Persistence of unvisited sites in presence of a quantum random walker 
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A study of persistence dynamics is made for the first time in a quantum system by considering the 
dynamics of a quantum random walk. For a discrete walk on a line starting at x = at time t = 0, 
the persistence probability P(x, t) that a site at x has not been visited till time t has been calculated. 
P(x,t) behaves as (t/\x\ — l)~ a with a ~ 0.3 while the global fraction V(t) = P(x,t)/2t of sites 
remaining unvisited at time t attains a constant value. F(x,t), the probability that the site at x 
is visited for the first time at t behaves as (t/\x\ — I) -13 /\x\ where /3 = 1 + a for t/\x\ » 1, and 
T(t) — 53 F(x,t)/2t ~ A few other properties related to the persistence and first passage 

times are studied and some fundamental differences between the classical and the quantum cases 
are observed. 

PACS numbers: 05.40.Fb,03.67.Hk 



Various studies on persistence in dynamical systems 
have been made in recent years [l[ where the dynami- 
cal process is purely classical in nature. The persistence 
probability that the order parameter in a magnetic sys- 
tem has not changed sign till the present time [2| and 
the persistence of unvisited sites in a diffusion problem 
[H are common examples which have received extensive 
interest. The importance of persistence phenomena lies 
in the fact that the persistence probability in many sys- 
tems shows an algebraic decay with a novel exponent. 

The dynamics of a quantum system is expected to be 
different from the corresponding classical case. However, 
to the best of our knowledge no attempt to study persis- 
tence in quantum systems have been made so far. Here 
we report a study of persistence and related dynamical 
features in a quantum system which has a classial ana- 
logue such that it is possible to determine the differences, 
if any, between them. 

The dynamical process we consider is a discrete quan- 
tum random walk (QRW). Classical random walk (CRW) 
on a line is a much studied topic 0, [H, @ where at ev- 
ery step one tosses a fair coin and takes a step, either 
to the left or right. In a quantum generalisation of a 
classical walk, one may consider superposition of particle 
movement in both directions, however such a process is 
physically impossible due to non-unitarity. A QRW in 
one dimension thus involves a quantum particle on a line 
with an additional degree of freedom which may be called 
chirality @, Q ■ The chirality takes values left and right 
analogous to Ising spin states ±1, and directs the motion 
of the particle. At any given time, the particle is in a 
superposition of left and right chirality. Thus the wave 
function has two components and at each time step, the 
chirality undergoes a rotation (a unitary transformation 
in general) and the particle moves according to its final 
chirality state. 

The reason to study the qunatum random walk is that 
it is already known to have some distinctive behaviour 
when compared to the classical walk. While in the CRW, 



after t time steps, the variance a 2 oc t, the quantum 
walker propoagates much faster with a cx t. This feature 
of faster flow of a quantum walker has made its applica- 
tion in quantum computation highly relevant. Moreover, 
the QRW spreads roughly uniformly over the interval 
[— t/v2) t/y/2] whereas in the classical case the distribu- 
tion is peaked at the origin and drops off exponentially 
fast. 

In the case of the classical random walker, the well 
studied quantities are (a) the first passage time F c i(x,t) 
which is the probability that the walker has reached the 
site x for the first time at time t [6( and (b) persistence 
or the survival probability P c i(x, t) defined as the prob- 
ability that the site at x has not been visited till time t 
d, H . A related quantity is the average persistence prob- 
ability given by J2 X Pci{x, t)/L (where x can take values 
from — L/2 to L/2). This quantity is global and can be 
interpreted as the fraction of sites remaining unvisited 
till time t. 

In this paper, we have studied both the first passage 
probabilities and the persistence probabilities of sites in 
presence of a quantum random walker. Not only do we 
find out F (x, t) and P{x, i) as functions of t and x, we 
also observe several other interesting features of the dis- 
tributions which indicates that there are two independent 
dynamical exponents for the quantum random walker. 

In this context, it may be mentioned that in the QRW, 
the walks can be realized in several possible ways, e.g., 
a walk can be infinite timed in which the walk is al- 
lowed upto time t at which all measurements are made. 
Another method may be to have a semi infinite walk 
where the question is asked at every time step whether 
the walker is in a specific location, and if it is so, the 
system is allowed to collapse to that state. In our case, 
we simply evaluate the occupation probabilities at each 
step for each discrete location and calculate the persis- 
tence and first passage probabilities which are related to 
it. There is no actual measurements being made which 
will lead to a collapse to any particular state. Rather, 
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our calculations will correspond to what would be the 
result of repeated experiments. 

The states for a quantum random walker are written 
as |a;)|d) where x is the location in real space and the d 
the chirality having either "left" and "right" values. The 
chiral states are denoted as \L) and \R). A conventional 
choice of the unitary operator causing the rotation of the 
chirality state is the Hadamard coin represented by 



(i) 



The rotation is followed by a translation represented 
by the operator T: 



T\x)\L) 
T\x)\R) 



\x-l)\L) 
\x+l)\R) 



(2) 



The two component wave function ip(x,t) describing 
the position of the particle is written as 



ip L {x,t) 



and the occupation probability is given by 

f{x,t) = \i> L {x,t)\ 2 + \i> R {x,t)\ 2 . 



(3) 



(4) 



Normalisation of the probability implies that at every 
ste P Ex/OM) = !• 

We have used two methods to evaluate f(x,i): 

1. A quantum random walk is generated using the above 
operators and ipL(x,t) and ipn(x, t) are evaluated numer- 
ically for all x and t. We start with the walker at the 
origin with V>l(0, 0) = ao,ipn(0,0) = bo; Oq + &q = 1 and 
all other ipL and tpR taken equal to zero. 

2. In the second method we use the expressions Q 



i + (-ir +t rdk 

1pL{X,t) = / — (1 



cos k 



2tt Vl + cos 2 k 



-i(ujkt+kx) 



(5) 



-l) x+t f dk 



-i(uikt+kx) 



27r vl + cos 2 k 



(6) 

(which are obtained for a initial state with left chirality, 
i.e., ao = 1, bo = 0) and evaluate f(x,t) directly. 

The persistence probability that a site x has not been 
visited till time t, is given by in terms of the occupation 
probabilities: 



p( x ,t)= n (i-f(x,t)), 



while the first passage time 



F(x,t)= -f(x,t'))f(x,t) 



(7) 



(8) 
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FIG. 1: The data collapse for the persistence probability and 
scaled first passage probability for a quantum random walk 
shown for different values of x. 
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FIG. 2: Plot of average persistence V(t) and first passage 
time T{t) indicate V(t) is approaching a constant for large 
time whereas T{t) behaves as t~ l . 



= P(x,t-l)f(x,t). 

Some time dependent features like hitting time, recur- 
rence time etc. @, [T3| of the QRW have been studied 
earlier which involve the first passage time. However, in 
these studies, the spatial dependence has not been con- 
sidered. For example, quantities like first passage time 
specifically at the origin has been estimated. 

In the first method of calculation, in general, real val- 
ues of the coefficients ao and bo give an asymmetric re- 
sult for f(x,t) in the sense the occupations probabili- 
ties f(x, t) ^ f(—x, t) as there is destructive interference. 
However, taking specific complex values of ao and 6o or 
by taking a different coin results in a symmetric walk 0] . 

We first discuss the results for the symmetric case ob- 
tained by taking a = l/V2,&o = i/V%- We find that 
P(x,t) and F(x,t) as functions oft and x behave as 

P(x,t)^(t/\x\^l)- a (9) 

for t/\x\ > 1, and 

F{x,t)<x(t/\x\-i)-P/\x\ (10) 

for t/\x\ » 1 with a ~ 0.31 and (3 ~ 1.31. The data 
collapse for P(x, t) and F(x, t) are shown in Fig. [TJ 

Using the above form of P(x,t), the average proba- 
bility V{t) that a site is unvisited till time t is given 
by ^ _. P(x,t)/2t (since x can assume discrete values 
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FIG. 3: Typical variation of F(x,t) against x for different 
values of t. The peaks are approximately fitted to 0.8lx m ° a ^ 9 . 



FIG. 4: Typical variation of F(x, t) against t for different 
values of x. The peaks are approximately fitted to 
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from —t to +t) which in the symmetric case is approxi- 
mately equal to 

- f P(x,t)dx = ~ f (--l)- a dx. (11) 
t Jo ' Jo x 

The integral on the rhs of eq (fTTj) can be done exactly 
giving a constant value equal to 0.738 (using A = 0.626 
which we get from the fitting of P(x, t) shown in FigfT]), 
indicating that V[t) is time independent. We indeed 
find by numerically evaluating ^2 X=0 P(x,t)/t that it ap- 
proaches a constant value as t becomes large. The con- 
stant is w 0.66, quite close to the value obtained above. 
On the other hand, Y2x=-t F(x, t)/2t varies as t~ 1 , a re- 
sult which can be obtained using a similar approximation 
and also obtained numercially. The numerical results are 
shown in Fig. [2] 

The function F(x,t) in fact has some additional fea- 
tures. Plotting F(x, t) against x or t, we notice that it 
has an oscillatory behaviour. These oscillations which 
die down for large values of t/x as is apparant from 
Fig. Q] can be traced to the oscillatory behaviour of 
f(x,t) for a QRW observed earlier @, H]. From Figs [3] 
and [U we observe that F(x,t) actually attains a max- 
imum value F max (x,t) at values of \x\ = x max (or 
t = t max ) for fixed values of t (or x). We notice that 
t) oc x max where 6 ~ 0.59. Keeping x fixed, 
F ma x{x, t max ) versus t max shows the same kind of depen- 
dence, i.e., F max (x, t max ) oc t m s ax . That the scalings with 
tmax and x m ax turn out to be identical is not surprising 
as x scales as t in a QRW. It is not possible to obtain this 
scaling form directly from eq. (fTU)) since F(x, t) attains 
a maximum value when t/\x\ is close to unity where the 
fitted scaling form is not exactly valid. In fact, eq. ([TO)) 
does not give any maximum value at all. 

Another dynamic quantity called hitting time has been 
estimated earlier for the QRW, in which an absorber is 
assumed to be located at a specific vertex of a hypercube 
within which the walk is conceived Q. The average hit- 
ting time is by definition the average time to reach that 
particular vertex for the first time. One can evaluate 
the average hitting time 77,(0;) using 77, (x) = Y^o tF(x, t) 



where t is allowed to vary from to T: 

r h w f tF(x,t)dt 
Jo 

- T 2 - f3 x (} - 1 /(2 -0) + 0{T- p+1 ). 

The numerical data (not shown) gives a fairly good agree- 
ment with this scaling. The above equation shows that 
Th blows up for T — > 00 in agreement with some earlier 
results using other coins [§| . 

When ao and bo are allowed to take up real values, in 
general the probability f(x, t) is asymmetric. Now we 
study the behaviour of P(x, t) and F(x, t) for x > and 
x < independently. In particular, we find that when 
there is a left bias, i.e., the probability that /(x < 0, t) > 
f(x > 0,i), P(x < 0,t) is less than P(x > 0,t), which 
is to be expected. Moreover, the scaling behaviour for 
x < shows an additional dependence on x, P(x < 0,i) 
is now given by 

P{x,t) oc (t/\x\ - l)- a \x\ p 

where p is approximately 0.03 for ao = 1, bo = 0. If a and 
b are chosen differently, the value of p shows up a vari- 
ation with the chosen values, however it is still 0(0.01). 
Such a small value of p suggests that this could be due to 
numerical errors. Fig. 5 shows the left-right persistence 
behaviour. 

The first passage probability for x < shows a similar 
correction to scaling 

F(x,t) oc (t/\x\ - l)-P\x\- q 

where q = 1 — 0(0.01) Fig. 6 shows typical variations of 
F(x, t) for a = 1, b = with q = 0.93 . 

The significance of the various results obtained in the 
present work becomes quite clear when these are com- 
pared to those of the classical random walker. For a 
CRW, the first passage time F c i(x,i) is known exactly 

\x\ 

Fd{x,t) = const -^exp(-x 2 /£)t), 
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FIG. 5: Data collapse of P(x,t) for a quantum random walk 
with a left bias. For the right side data, p — and for the left 
side, p — 0.03. 
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FIG. 6: Data collapse of F(x,t) for a quantum random walk 
with a left bias. For the right side data, q — 1 and for the left 
side, q = 1 - 0.07 



where D is a diffusion constant. Hence F(x, t) behaves as 
i~ 3 / 2 for large t. The persistence or survival probability 
is given by erf (^J^^) which behaves as | | i 1 / 2 as t — ► oo. 
In fact the two quantities, persistence probability and 
first passage time are related in the continuum limit: 



F(x,t) = - 



dP(x,t) 

at ' 



(12) 



and hence the scaling behaviour (with time) of F can be 
found out from that of P or vice versa; e.g., (3 c i = ay + 1, 
where P c i cx t~ a " 1 and F c i cx t" 13 " 1 . 

In general one can write x oc t 1 for these kind of dy- 
namical phenomena where the inverse of 7 is a dynamical 
exponent. It can be expected that the scaling behaviour 
of all other quantities will be dictated by the exponent 7. 



Here we would like to emphasise that while this is true 
for the classical case, for the QRW, that is not the case. 
Let us consider each physical measure for the CRW and 
QRW to establish this claim. 

In the CRW, 7 = 7 C ; = 1/2 and the scaling of the first 
passage time with time occurs with an exponent j3 c i = 
37c/ and the persistence probability with an exponent 
&d = Jch For the quantum walker, 7 = j q = 1, however, 
the persistence behaviour and the first passage times vary 
in time with exponents a and /? which are not simple 
multiples of -f q although F(x, t) and P(x, t) indeed obey 
the relation given in eq. (p~2|) such that (3 = 1 + a. The 
form of the persistence probability and first passage time 
are in fact entirely different from those of the classical 
case as indicated by the collapsed data. 

In the classical case, the global fraction T c i behaves 
as t~ x . From this, one can infer that the global fraction 
Vd is approximately a constant. In the quantum case, 
surprisingly, we have identical scaling behaviour for P(t) 
(which saturates to a constant value) and T(t) (which 
varies as t~ r ). In terms of 7 however, T c i oc i~ 27e! while 
T cx t~ lq . In the QRW a striking feature is that these 
global properties do not involve the exponents a and f3 
at all and only j q dictates their behaviour. 

Concerning the maximum values of the classical prob- 
ability F c i , we find that it behaves as 1 /t max (for x con- 
stant) or ^jx 2 max (for t constant) showing that the ob- 
tained exponents are simple multiples of 7 C ; = 1/2. On 
the other hand, the behaviour of F max in the quantum 
case appears to depend on the value of a and not j q as 
it varies with t max or x max with an exponent 8 which is 
very close to 2a numerically. 

The average hitting time for a CRW is found to vary as 
T 7cl . In the QRW, this variation is given by T 2-/3 . For 
the classical case, 2 — (3 c i = 7 C ; but since no such relation 
exists for the quantum case, the hitting time scaling is 
therefore not dictated by "f q but by f3 (or a) only. 

Thus the scaling forms of different quantities for the 
quantum and classical walks are not identical in general 
and the scaling behaviour in the quantum case is dictated 
either by 7 9 or by a. From this we draw the important 
conclusion that there are two nontrivial and independent 
exponents for the QRW in contrast to the CRW. This re- 
sult may have serious impact on quantum computational 
aspects. 
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